cholesky_rank1_update Interface

public interface cholesky_rank1_update

Module Procedures

private pure subroutine cholesky_rank1_update_dbl(r, u)

Computes the rank 1 update to a Cholesky factored matrix such that .

Arguments

Type IntentOptional Attributes Name
real(kind=real64), intent(inout), dimension(:,:) :: r

On input, the N-by-N upper triangular matrix . On output, the updated matrix .

real(kind=real64), intent(inout), dimension(:) :: u

On input, the N-element vector . On output, the rotation sines used to transform to .

private pure subroutine cholesky_rank1_update_cmplx(r, u)

Computes the rank 1 update to a Cholesky factored matrix such that .

Arguments

Type IntentOptional Attributes Name
complex(kind=real64), intent(inout), dimension(:,:) :: r

On input, the N-by-N upper triangular matrix . On output, the updated matrix .

complex(kind=real64), intent(inout), dimension(:) :: u

On input, the N-element vector . On output, the rotation sines used to transform to .